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1. Introduction 


The origin and evolution of continents represent one of the greatest challenges in the field of Earth Science 
research today. Unraveling the history, internal dynamics, and evolution of the lithosphere, both in continental 
and oceanic settings, can provide insights into how the Earth's crust formed, the inception of plate tectonics, and 
the development of habitable continents. In the modern Earth, this knowledge enhances our understanding of the 
nature and evolution of tectonic plates [1], The present-day configuration of the Western passive continental 
margin of India is a consequence of its extensive and intricate tectonic history. 


The Western Ghats, known as the Sahyadri in its present-day configuration, ranks among the world's most 
spectacular and extensive escarpments, alongside those found in Australia, Brazil, and the South African passive 
margins. This long-uplifted topographical scarp (Fig. 1b) represents a prominent geomorphic feature. It 
encompasses uneven topography, deep valleys, numerous impressive waterfalls, dense forests, a complex 
geological history, extensive volcanic activity, and various structures and petro-physical properties [2]. It stretches 
approximately 1600 kilometers in a nearly straight line, running in a north-northwest to south-southeast direction. 
The study area, with an average elevation of 1200 meters, traverses across six states: i) Gujarat ii) Maharashtra iii) 
Goa iv) Karnataka v) Kerala, and vi) Tamil Nadu. The maximum elevation reaches 2500 meters, and the 
escarpment is aligned roughly parallel to the west coastal region (Arabian Sea) of Peninsular India. Along the 
Western Ghats, there is only one significant low-elevation break known as the Palghat Gap. It spans approximately 
30 kilometers in width, extending from the northern region of Gujarat to Kanyakumari. 


The evolution of the topographic scarp in the Western Ghats has been studied extensively, with recognized 
perspectives from various fields including geomorphological, geochemical, geological, and geophysical studies. 
The most popular explanations are marginal rift flank uplift [3-5], mantle plume-related drifting [6, 7], magmatic 
under-plating [8, 9], lithospheric delamination [10], isostatic compensation [11, 12], increasing tectonic processes 
and denudation rates [13-18], and flexural adjustment [19-21]. As per [22], the causes of upliftment can be 
categorized into three distinct ways: (i) Density reduction, either mechanically or thermal processes, in response to 
isostatic adjustments. (ii) Crustal buoyancy, attributed to an increase in lithospheric thickness, and (iii) Plastic 
necking, resulting from stretching and uneven denudation of lithosphere on both sides of the scarp. 


Recent studies have put forth various hypotheses regarding the formation of the Western Ghats. These include 
active rifting and extensional tectonics [23], igneous underplating [24], flexural isostatic adjustment of uplifted 
topography, later modified by tectonic and denudation processes [25], crust and mantle interactions during 
tectonic-magmatism associated with the Deccan volcanic provinces [26], and multistage rifting episodes [27-29]. 


While numerous studies, as mentioned earlier, have proposed multiple hypotheses and theories to 
comprehend the evolution of the WG and its connection with the expansion of the Western Passive Continental 
Margin (WPCM), a conclusive explanation for the evolution of the WG remains elusive [30, 31]. This is due to the 
region's diverse lithological variations, complex interactions between the crust and mantle, the presence of the 
Palghat gap interrupting the topography, and the exposure of lower crustal rocks in the entire Southern Granulite 
Terrain (SGT) [32]. Here, we conducted a comparative analysis of various models of crustal and lithospheric 
thickness in the study region. We aim to evaluate the effectiveness of different modeling techniques and the 
influence of varying data sources on these models' quality. Additionally, we have developed 2-D lithospheric 
density models along the Western Ghats (WG) using satellite-derived EGM 2008 gravity data. These models have 
been constrained by previous results, aiming to enhance our understanding of the factors contributing to the 
evolution of the WG by enabling us to decipher the underlying crust and lithospheric structure. 


2. Geology and Tectonic Setting 


The Western Ghats are a highly topographic region in India, after the Himalayan Mountain building. The 
present WPCM of India and the surrounding regions formed through several stages of rifting episodes that 
occurred before the break-up Gondwana supercontinent around 180 Ma [33, 34]. Geological history indicates that 
during the early Jurassic period, the Gondwana supercontinent started to divide into two major landmasses: 1) 
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West Gondwana and 2) East Gondwana, as a result of continental drift. The West Gondwana is comprised of South 
America and Africa, which consists of a mosaic of Archean/Paleoproterozoic cratons and Neo-Proterozoic to 
Cambrian Pan-African orogenic systems [34-36]. On the other hand, the East Gondwanaland drifted south-south 
easterly relative to Africa in the Late Jurassic (170-120 Ma) constituted by India, Madagascar, Seychelles, Sri Lanka, 
Australia, and East Antarctica. Furthermore, within the West Gondwana, India and Madagascar began drifting 
away from Antarctica-Australia at approximately 140 Ma [37], and the dispersal process continued. [38] suggested 
that this could be attributed to thermo-mechanical interactions between these two continental lithospheric 
fragments. Later, at approximately 90 Ma, the India-Seychelles landmasses began to rift rapidly, moving 
northward and eventually separated from Madagascar due to the activity of the Marion plume [39] (Reeves, 2014 
and references therein). At approximately 65 Ma, the Deccan flood basalts erupted as a result of the reunion 
plume hotspot activity, leading to the separation of the India-Seychelles continental blocks [40]. In general, 
continental lithospheric plate rifting and break-up are often associated with mantle plume activity. Consequently, 
many researchers believe that these two major rifting events strongly influenced the evolution and enlargement 
of the WPCM. 


Based on geological formations, the entire Western Ghats region (Fig. 1a) is divided into three parts: Northern, 
Middle, and Southern. The northern part, extending from latitude 16° to 20° N, is predominantly occupied by the 
Deccan Volcanic Provinces (DVP) and is commonly referred to as the Konkan coast. The Deccan Traps region is 
primarily characterized by thick sequences of volcanic lava flows, approximately 2 kilometers in thickness, 
concealing Mesozoic sediments and the remnants of Deccan volcanism beneath this zone [41]. The middle zone, 
situated between latitudes 12° to 16° N, is characterized by Precambrian volcano-sedimentary sequences. The 
Dharwar Craton in Peninsular India is one of the oldest continental Archean crusts, primarily composed of 
tonalite-trondhjemite gneisses and granites [42]. These rocks range from ~3600-2500 Ma and overlie the 
basement complex gneisses. Furthermore, the Dharwar Craton has been subdivided into two tectonic blocks: the 
Western Dharwar Craton (WDC) and the Eastern Dharwar Craton, based on differences in age, deformation, and 
tectonic history. The southern part, known as the Malabar Coast (latitude extending from 8° to 12° N), consists of 
high-grade granulite rocks found within the Southern Granulite Terrain (SGT). This highly metamorphosed SGT 
represents one of the world's most deeply exposed Precambrian continental crusts, which underwent 
deformation approximately 3500-550 Ma [43]. 


3. Crust and Lithospheric Variations 


Understanding crust and lithospheric variations is crucial for comprehending the internal structure and 
geodynamic evolution of tectonic plates, as well as the formation of elevated surface topography. Global studies 
on crust and lithospheric thickness have revealed significant variability like the Earth's, comprising the crust and 
upper mantle, across different tectonic and lithological regions. In some locations, this boundary is well-defined, 
while in most others, it exhibits a diffused nature [44]. Crust and lithospheric thickness have been estimated using 
various geophysical methods over the past two decades in the Western Ghats and surrounding areas. In this 
region, several studies, either in part or as a whole, have examined the estimated crust and lithospheric thickness. 
These investigations have been conducted by various researchers using different geophysical methods, which are 
explained below. 


3.1. Crust 


From gravity studies, crustal thickness following range: SGT (35-44 km), EDC (32-38 km), WDC (40-52 km), DVP 
(30-50 km), and along the WPCM (38-53 km); [25, 45-49]. Various seismic and receiver function analyses: SGT (36- 
46 km), EDC (31-38 km), WDC (40-60 km), DVP (30-54 km), and along the WPCM (37-45 km); [26, 29, 50-56]. 


3.2. LAB 


The estimated LAB depths by different techniques (Surface wave analysis, receiver functions, Magnetotelluric 
data, gravity, and geoid anomaly) made by various authors is as follows: SGT (120-185 km), EDC (70-185 km), WDC 
(100-200 km), DVP (95-120 km), and beneath the WPCM (~160 km); [57-64]. 
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4. Data Utilization 


Geophysical potential field methods, particularly the gravity method, measure the Earth's gravitational field as 
influenced by surface and subsurface rocks, structures, and physical properties such as density distributions. 
Satellite-derived global gravity models have limitations because they do not represent short wavelengths in the 
Earth's gravity field. In this study, we have utilized various datasets to comprehend subsurface density variations 
and provide insights into the geodynamic prospects beneath the Western Ghats. 


4.1. Topographic Map 


The WPCM of India forms an uninterrupted coastal wall that spans from 8° to 22° N in latitude. It stretches 
over a length of approximately 1600 km, extending from the Tapti River in the north to Kanyakumari in the south. 
The only interruption in this expanse occurs in the Palakkad region near the southernmost part of India, where it 
narrows to roughly 30 km in width. This study used from the Shuttle Radar Topography Mission (SRTM) 15+ for 
topographic analysis [65]. The topography map (Fig. 1b) depicts all major geological formations and exhibits a 
strong correlation with the regional geology (Fig. 1a) in the study area. In the coastal region gradually ascends ina 
step-like pattern, marked by significant elevation changes. The topographic variations range from a minimum of - 
3.70 km to a maximum of 2.60 km. High topographic values are predominantly concentrated along the 
escarpment, while the oceanic region exhibits the lowest values. In particular, the southern part of the WDC and 
the northern section of the SGT feature prominently elevated topography, with the maximum elevation reaching 
2.6 km in the Nilgiri block. In contrast, the EDC exhibits a flatter topography, averaging around 0.3 km, in 
comparison to other regions such as the DVP, WDC, CB, and SGT. 


4.2. Complete Bouguer Anomaly Map 


In the present study region, research efforts have been limited due to the challenging terrain characterized by 
hilly landscapes and a lack of accessible roads. This has posed significant logistical challenges for conducting on- 
site investigations. As a result, we have utilized satellite-based measurements derived from the Earth Gravitational 
Model (2008) to overcome these limitations and gain insights into crustal and lithospheric variations across the 
Western Ghats. The EGM2008 model [66] is a spherical harmonic model, which is obtained from satellite missions, 
including CHAMP and GRACE [66]. This model is freely accessible and contains coefficients extended up to a 
degree of 2190 with an order of 2159. The gravity method measures changes in the earth's gravitational field on 
the earth's surface to determine the subsurface lateral density variations, providing valuable indirect information 
on the mass distribution at depths. 


The Complete Bouguer Anomaly (CBA) map, which has been terrain corrected to account for the influence of 
topography, is shown in Fig. (1¢). This map provides a visual representation of gravity values ranging from -124 to 
168 mGal. Gravity anomalies, both positive and negative, are depicted across the study area, offering insights into 
the subsurface density variations and mass distribution beneath the Earth's surface. Positive anomalies are 
predominantly observed in the oceanic and coastal plain regions, indicating areas with relatively higher mass 
concentrations beneath the Earth's surface. Conversely, negative anomalies are more prevalent in continental 
regions, with the largest negative anomalies observed in the western continental parts of the DVP, WDC, SGT, and 
CB. These negative anomalies are often associated with the presence of thick crustal roots. Remarkably, over the 
WPCM, we observe a distinct negative anomaly, especially in the NNW-SSE direction within the WDC. This region 
exhibits the lowest gravity anomalies, ranging from -124 to -35 mGal. Notably, this anomaly pattern closely 
correlates with the topography map (Fig. 1b), where the areas of low gravity anomalies align with the elevated 
topographic features along the WPCM. 


5. Radially Averaged Power Spectral Analyses 


Radially Averaged Power Spectral Analysis is a commonly used statistical method for interpreting gravity and 
magnetic anomalies [67] to calculate the average depth of interface layers caused by subsurface sources. In 
this study, we applied 2-D Fourier Transformation analysis to the Complete Bouguer Anomaly (CBA) to distinguish 
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Figure 1: a) Regional geology map of the study area (after, Geological Survey of India, 1995). SGT: Southern granulite terrane, 
WDC: Western Dharwar Craton, CB: Cuddapah Basin, EDC: Eastern Dharwar Craton, EGMB: Eastern Ghat Mobile Belt, GG: 
Godavari Graben, DVP: Deccan Volcanic Province, BC: Bastar Craton, CG: Closepet Granite. Geology and tectonic map of study 
area shown grey coloured in rectangular box. b) Topography map of the study area (in km). The dashed yellow lines indicate 
previous gravity studies from various authors, the blue solid line represent escarpment. The references with red-colored 
backgrounds labeled as (A, B, C, D) have been selected to estimate the uncertainties in crust and lithospheric thickness. c) 
Complete Bouguer anomaly (in mGal) map. 


between long-wavelength and short-wavelength signals associated with deeper and shallower sources, 
respectively. The CBA power spectrum plot, with wavenumber on the x-axis and the logarithm of spectral energy 
on the y-axis, shows a continuous decrease with increasing wavenumber. The gradient of the linear segments in 
the Fourier power spectrum is related to the average depth of various subsurface layers. Using the slopes of these 
lines, we calculated mean depths for the interface boundaries, which are denoted as L1= 176 km (LAB), L2= 47 km 
(Moho), L3= 21 km (Middle crust), and L4= 10 km (Upper crust) (Fig. 2). 
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Figure 2: Radially averaged power spectrum calculated averaged depth to sources Complete Bouguer anomaly estimated 

different layers depths shown red colours. L1= 176 km (LAB), L2= 47 km (Moho), L3= 21 km (Middle crust), and L4= 10 km 

(Upper crust). 
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6. Uncertainties in Crust and Lithospheric Thickness 


Many geophysical methods are non-unique, as they utilize various physical parameters to gain insights into 
and investigate the Earth's subsurface. The challenge lies in producing reliable geological images, primarily due to 
the limited data sources available on the Earth's surface, including deep seismic sounding, potential field data, 
receiver functions, and seismic tomography. Consequently, uncertainties arise during the interpretation of 
geophysical data, making it challenging to reach definitive conclusions. In this study, we interpreted a total of 12 
profiles based on previous research results, incorporating a diverse set of gravity and magnetic anomalies, along 
with global models. Our objective was to estimate uncertainties in both crustal and lithospheric thickness. We 
employed the following formula to calculate these uncertainties: 


Uncertainty formula = V[Xi(x; — 1)2/(n x (n — 1))] 
In this formula, x; represents individual data points, u is the mean, and n is the number of data points. 


e The West-East-oriented profile (Fig. 3a) along the 9.5°N latitude, which spans approximately 1000 km across 
the Arabian Sea, SGT, and the Bay of Bengal, clearly reveals positive gravity anomalies over the ocean and 
negative anomalies in the continental region. While the geoid anomaly does not exhibit a strong correlation 
in the western ocean part, it does align well with elevated areas. Regional gravity and isostatic anomalies 
indicate positive signatures, suggesting a shallower crustal and lithospheric depth, whereas negative 
anomalies imply thicker structures. According to [68], who used an integrated modeling approach, the 
Moho (crust-mantle boundary) and LAB (Lithosphere-Asthenosphere Boundary) depths are estimated to be 
around 12-15 km over the ocean and 27-36 km beneath the SGT. This depth range closely matches the 
findings of [69] in their global model (75%), but it deviates significantly from other global models and 
receiver function results, which display substantial variation. The LAB depth beneath the ocean ranges from 
67 to 128 km, while it's approximately 130 km under the SGT, and these values do not align with any 
existing models. The maximum observed uncertainties in crustal and lithospheric thickness along this 
profile are 7.75 km and 18.79 km, respectively, with average estimated uncertainties in the crust and 
lithospheric variations at 13.04% and 23.92%. 


e In another West-East profile (Fig. 3b) along the 13.5°N latitude, the gravity anomaly exhibits a continuous 
decrease from west to east, except at the end of the profile, closely mirroring the variations in Moho and 
LAB depths. The Moho values obtained from [26] align closely with the trends observed in two global 
models, Crust1.0 and [69], but only match more than 85% with receiver function results over the WDC, CG, 
and EDC regions. The LAB depths observed beneath the WDC and EGMB indicate thickness, while those 
beneath CG and EDC appear thinner, deviating from the patterns seen in global models. The maximum 
observed uncertainties in crustal and lithospheric thickness along this profile are 5.86 km and 18.67 km, 
respectively, with average estimated uncertainties in the crust and lithospheric variations at 20.01% and 
21.80%. 


e Vasanthi and Santosh [70] estimated lithospheric structure based on satellite-derived gravity data (Fig. 3c). 
Their results, obtained using the finite element technique, revealed clear demarcations of the boundaries 
of major geotectonic formations in the Dharwar craton, including the western, central, and eastern regions, 
as well as the surrounding areas. The averaged depths of the Moho and LAB beneath the various major 
blocks within the Dharwar craton were as follows: 38 km and 138-175 km for the western block, 41 km and 
102-138 km for the central block, and 42 km and 88-120 km for the eastern block, respectively. At the 
boundaries between the western and central blocks, the depths were approximately 52 km and 147 km, 
while between the central and eastern blocks, they were around 46 km and 103 km. Notably, the results of 
this study did not align with any global models or receiver function analysis. The maximum observed 
uncertainties in the crust and lithospheric thickness along this profile were 7.85 km and 21.26 km, 
respectively. The average estimated uncertainties in the crust and variations in lithospheric thickness were 
7.82% and 12.85%, respectively. 


e A 1300 km-long profile in the SSW-NNE direction (Fig. 3d) traverses the SGT, WDC, and Cuddapah basin, 
revealing numerous gravity and magnetic anomalies, including highs and lows. The crustal thickness, as 
estimated by [71] Kumar et a/. (2020) using 2-D lithospheric density modeling, shows variations. It measures 
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Figure 3: Uncertainty analysis estimated maps from Previous studies correlation with available global and receiver functions 
results. A) Profile no.1 showing the yellow dot lines in (Fig. 1b) taken from Kumar et al., 2013. B) Profile no. 4 showing the 
yellow dot lines in (Fig. 1b) taken from Kumar et al., 2013. C) Profile no.5 showing the yellow dot lines in (Fig. 1b) taken from 
Vasanthi and Santosh, 2021. D) Profile no.11 showing the yellow dot lines in (Fig. 1b) taken from Kumar et al., 2019. SGT: 
Southern granulite terrane, WDC: Western Dharwar Craton, CB: Cuddapah Basin, EDC: Eastern Dharwar Craton, EGMB: Eastern 
Ghat Mobile Belt, GG: Godavari Graben, CG: Closepet Granite. 


approximately 25 km in the ocean, around 41 km under the SGT, and between 38-40 km in the WDC and 
Cuddapah basin regions. The calculated lithospheric thickness exhibits variations beneath the SGT and 
Eastern Indian shield, measuring between 120-130 km. In contrast, the WDC region displays a lithospheric 
thickness of 160-180 km. These values do not correlate with global models or receiver function data. The 
maximum observed uncertainties in the crust and lithospheric thickness along this profile are 4.27 km and 
24.18 km, respectively. The average estimated uncertainties for crust and lithospheric variations are 18.47% 
and 14.26%, respectively. 


7. 2-D Lithospheric Density Modeling along the Western Ghat Escarpment 


Furthermore, we conducted 2-D lithospheric density modeling along the Western Ghat Escarpment (WPCM). 
This involved generating cross-sections along three lengthy profiles, each spanning approximately 1600 km and 
oriented in the SSE-NNW direction. These profiles, labeled as A-A’, B-B', and C-C' (Fig. 4a), traverse various 
lithological blocks within the region, including the SGT, WDC, and DVPTo perform this modeling, we utilized 
Geosoft GM-SYS software, originally developed for potential fields by [72, 73]. It enabled us to perform forward 
modeling effectively. In a Cartesian coordinate system, the gravity anomaly profile of a 2-D anomalous density 
source in a plane perpendicular to the strike (in the y-directional extension) can be approximated using a polygon 
[72]. Improving the accuracy of this geometric representation involves increasing the number of vertices in the 
polygon. According to [72], The vertical component of the gravity field anomaly (g,) at point P (xi, zi) due to an 
infinitely extending mass source with a density contrast (Ap) in the y-direction can be expressed as: 


_ Ap. (z — 2) 
92% Zj) = 2G ff @=ay eG oo) dx dz 


Where G is the universal gravitational constant. 


It's important to note that interpreting potential data modeling often encounters non-uniqueness challenges. 
This is primarily due to the lack of direct information related to subsurface structures, petrophysical 
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characteristics, and the composition of various crustal and lithospheric layers. Given this complexity, it's possible 
to create an infinite number of models that produce the same observed gravity anomaly. To address this 
ambiguity, we employ a strategy of constraining the model using previously acquired complementary information. 


In our approach, we have incorporated available seismic constraints obtained from receiver functions in 
various studies to determine both crustal thickness and the depth of the LAB. These constraints have been drawn 
from a range of sources, including gravity studies with seismic constraints [74], converted wave techniques such 
as P and S receiver functions [60], 2-D lithospheric density modeling developed by [68], a 3-D lithospheric model 
established using surface wave observations and classical tomographic methods [61], as well as a finite element 
approach applied by [70]. The initial depths of the model's layers have been derived from radially averaged power 
spectrum analysis (Fig. 2). 


7.1. Profile - A-A’ 


The western section of this profile has been specifically chosen along the coastal area, characterized by a range 
of gravity anomalies, with values fluctuating between -55 mGal and 41 mGal. The most notable positive gravity 
anomalies have been observed beneath the DVP, reaching approximately 40 mGal. The gravity model, constrained 
by seismic data and applied along this profile, provides insights into crustal thickness variations. Within this model, 
the SGT region displays a crustal thickness ranging from 33 km to 42 km, the WDC region ranges from 39 km to 42 
km, and the DVP region exhibits variability from 31 km to 43 km (Fig. 4b). Likewise, the LAB's variation is evident in 
the SGT region, spanning from 57 km to 83 km, in the WDC region, ranging from 76 km to 80 km, and in the DVP 
region, extending from 79 km to 94 km. 


7.2. Profile - B-B’ 


Profile B-B' has been precisely positioned over the escarpment, exhibiting the most significant negative 
Bouguer anomaly compared to the other profiles. The range of anomalies varies from a minimum of -118 mGal to 
a maximum of -37 mGal. Two noteworthy features stand out in this profile. First, near the geological boundary 
between SGT and WDC, there is a wavelength of approximately 200 km with a fluctuating anomaly of around -60 
mGal. Second, near the boundary between EDC and DVP, there is a wavelength of approximately 175 km with a 
fluctuating anomaly of about 35 mGal. Crustal thickness estimates have been determined based on constrained 
seismic results, revealing values for SGT ranging from 34 km to 45 km, WDC ranging from 37 km to 48 km, and 
DVP ranging from 34 km to 42 km. Additionally, LAB variations have been estimated within the Western Ghat, with 
constrained data available from select locations within the major geological blocks. These estimates indicate 
depths for SGT ranging from 87 km to 160 km, WDC ranging from 137 km to 166 km, and DVP ranging from 92 km 
to 137 km (Fig. 4c). Notably, a thick LAB, measuring between 150 km and 166 km, has been observed beneath the 
northern part of SGT and the southern part of WDC, corresponding to a highly elevated region. 


7.3. Profile - C-C’ 


This profile has been chosen on the eastern side of the uplifted escarpment, revealing numerous undulations 
in gravity anomalies. The topographic variations are largely linear, spanning approximately 750 meters, except at 
the boundaries between SGT and WDC. The gravity anomaly along the profile exhibits an overall range from -97 
mGal to +6 mGal, encompassing three major geological units. We have determined variations in the crustal 
structure beneath these units, with SGT ranging from 36 km to 47 km, WDC from 42 km to 46 km, and DVP from 
36 km to 42 km. Regarding the depth of the LAB, we observe depths of 114 km to 142 km in SGT, 115 km to 148 
km in WDC, and 82 km to 115 km in DVP. Notably, most of the SGT and the southern parts of WDC display similar 
crustal and lithospheric variations, approximately measuring 45 km and 140 km, respectively (Fig. 4d). 


8. Conclusion 


The accurate characterization of the subsurface of geological structures primarily relies on geoscientific input 
data. Nevertheless, uncertainties within this input data have the potential to yield unreliable geological structures, 
thereby affecting the reliability of geoscientific findings. In this study, we present several profiles depicting variations 
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Figure 4: 2D density lithospheric modelling a) Topography map of the study area, dashed orange lines indicate selected three 
(A-A, B-B’, C-C’) SSE to NNW profiles. The solid blue line is representing the Western ghat’s escarpment b) Profile A-A’ profile 
selected western part of the escarpment c) Profile B-B’ profile selected over the escarpment d) Profile C-C’ profile selected 
eastern side of the escarpment. The major geological boundaries shown by thick vertical black solid lines. D: Density, SGT: 
Southern granulite terrane, WDC: Western Dharwar Craton, EDC: Eastern Dharwar Craton, DVP: Deccan Volcanic Province. 


in crustal and lithospheric thickness that exhibit significant differences from one another. These variations in 
crustal and lithospheric thickness range from 8% to 20% and 13% to 24%, respectively. Such disparities can exert a 
profound influence on tectonic and geodynamic interpretations. Our research encompasses the Western Ghats 
and the surrounding regions of SGT, WDC, and DVP, characterized by distinct structural and physical properties. 
The precise mechanisms governing the evolution of the Western Ghats and its connection with the maturation of 
the Western Ghat Coastal Margin (WPCM) remain enigmatic. Therefore, we have undertaken 2-D lithospheric 
density modeling of crustal and lithospheric structures along the Western Ghat escarpment, utilizing three 
selected profiles. These profiles have been constrained by previous geophysical studies and globally available 
models. The density model unveils a five-layered configuration, providing insights into crustal thickness over SGT 
(ranging from 33 km to 47 km), WDC (ranging from 37 km to 48 km), and DVP (ranging from 31 km to 43 km). The 
lithospheric-asthenospheric boundary depths vary in SGT (ranging from 57 km to 160 km), WDC (ranging from 76 
km to 166 km), and DVP (ranging from 79 km to 137 km). The conspicuous changes in crustal and lithospheric 
structures along the Western Ghats escarpment hint at distinct tectonic and drifting episodes along the western 
margin of India. Further investigations are required to validate this concept and reduce uncertainties, contributing 
to a better understanding of the tectonic evolution of India's Passive Continental Margin. 
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